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We present a high-resolution photoemission study on the strongly correlated Ce-compounds 
CeCue, CeCu2Si2, CeRu2Si2, CeNi2Ge2, and CeSi2. Using a normalization procedure based on 
a division by the Fermi-Dirac distribution we get access to the spectral density of states up to an 
energy of 5fcsT above the Fermi energy Ef- Thus we can resolve the Kondo resonance and the 
crystal field (CF) fine-structure for different temperatures above and around the Kondo temperature 
Tk- The CF peaks are identified with multiple Kondo resonances within the multiorbital Anderson 
impurity model. Our theoretical 4/ spectra, calculated from an extended non-crossing approxi- 
mation (NCA), describe consistently the observed photoemission features and their temperature 
dependence. By fitting the NCA spectra to the experimental data and extrapolating the former to 
low temperatures, Tk can be extracted quantitatively. The resulting values for Tk and the crystal 
field energies are in excellent agreement with the results from bulk sensitive measurements, e.g. 
inelastic neutron scattering. 



PACS numbers: 71.27.+a 71.28.-|-d 79.60.-1 71.10.-w 
I. INTRODUCTION 

Since the experimental discovery of heavy-fermion 
(HF) compounds, metallic systems with rare-earth (4/) 
and actinide (5/) elements have been thoroughly investi- 
gated both experimentally and theoretically.^'^ The term 
heavy-fermions refers to the observation that these sys- 
tems behave as if the conduction electrons had an enor- 
mously high effective mass. The coefficient of the linear 
term in the low-temperature specific heat, the Sommer- 
feld coefficient, can be of the order of 1 JK~^mol~^, cor- 
responding to an effective electron mass enhancement of 
a factor 1000 compared to the free electron mass. In 
addition, these 4/ and 5/ systems show a variety of 
anomalous ground state and unusual low-temperature 
properties 

The single-impurity Anderson model^ (SIAM) embod- 
ies the key mechanisms of the many-body physics in 
heavy-fermion systems, namely a strong local Coulomb 
interaction of the /-electrons with the resulting local mo- 
ment formation and a weak hybridization between con- 
duction electron states and / states, which lead to the 



Kondo effect. ^'^ For cerium systems the theoretical de- 
scription is significantly simplified, because the Ce ions 
contain at most only one 4/-electron (i.e. 4/^). Em- 
ploying variational methods^'^, numerical renormaliza- 
tion group (NRG) calculationsi2iii or self-consistent dia- 
grammatic resummations, like the non-crossing approxi- 
mation (NCA) — and the conserving T- matrix approx- 
imation (CTMA)i^ii^ii^, one can calculate the / spectral 
density of states (DOS) near the Fermi level within the 
framework of the SIAM. The most important feature in 
the 4/ DOS of Ce-systems is a narrow peak with a max- 
imum just above the Fermi Energy Ep, the Kondo res- 
onance. The line width of the Kondo-resonance is given 
by the low-energy scale of the problem, defined by the 
Kondo-temperature Tk- For Ce systems kgTK is typi- 
cally a few meV, although the bare model parameters of 
the SIAM are usually two or three orders of magnitude 
larger. 

Although in the HF compounds the /-atoms form an 
ordered lattice, the simplest approach is the SIAM, em- 
ploying non-interacting local impurities, which has been 
very successful in even a quantitative description of high- 
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resolution photoemission spectral^. This can be under- 
stood as follows: The direct orbital overlap between ad- 
jacent /-atoms is negligible, so that the on-site Kondo 
resonances form coherent delocalized states below the lat- 
tice coherence temperature only due to their coupling via 
the conduction band. Therefore, their dispersion band 
width is essentially given by the Kondo resonance width 
itself, i.e. by the single- impurity Kondo temperature Tk- 
Since photoemission spectroscopy (PES) measures pre- 
dominantly the momentum integrated f-spectra, these 
PES spectra are well described by the SIAM. However, to 
describe the coherent heavy-fermion state with a dispers- 
ing narrow band close to the Fermi level, other models 
like the renormalized band picture^^ or the periodic An- 
derson model (PAM)i^ are required. 

PES, in particular angle resolved photoelectron spec- 
troscopy with excitation energies in the VUV-range 
(ARUPS), is one of the most direct experimental meth- 
ods to investigate the electronic structure of solids. For 
rare-earth systems, PES allows to study the occupied 
part of the 4/ spectral function with high accuracy. Un- 
fortunately, the main spectral weight of the Kondo res- 
onance appears above the Fermi level, where the photo- 
electron intensity is suppressed by the Fermi-Dirac distri- 
bution (FDD). On the other hand, inverse PES (IPES) 
can in principle measure the density of states above Ep, 
but the energy resolution — usually considerably larger 
than 100 meV — is not sufficient to investigate the details 
of the interesting 4/ spectral features, i.e. the Kondo 
resonance, the crystal field structures (CF), and in most 
cases even not the spin-orbit splitting. Recently it has 
been showniS, on CeCu2Si2 that PES with high energy 
resolution gives access to the Kondo resonance (KR). 
This first and direct observation of the KR was pos- 
sible by the application of a well known normalization 
procedure^ that allows to recover the thermally occu- 
pied DOS up to BksT above Ep. At higher tempera- 
tures this energy range can even cover the CF structures 
above the KR. From the normalized PES data one can 
determine the Kondo temperature Tk from the KR line 
width and the CF splittings of the J = 5/2 spin-orbit 
partner {Aqf)- A more detailed quantitative analysis 
is possible by comparing the experimental 4/ spectrum 
with theoretical spectral functions calculated by the non 
crossing approximation (NCA).^^ By iteratively fitting 
the data one can determine all relevant model parame- 
ters of the SIAM,^- from which eventually Tk and Acf 
can be extracted. 

In this paper we present a systematic investigation of 
the classical HE systems CeCug, CeCu2Si2, CeRu2Si2, 
CeNi2Ge2, and CeSi2, which show Kondo temperatures 
from 5 K to about 40 K as extracted from other ex- 
periments. As explained later, we restrict our investiga- 
tions to 7-Ce like materials showing a low hybridization 
strength and consequently low Kondo temperatures. We 
compare the PES results for the Kondo temperatures and 
the crystal field energies with values determined by bulk 
sensitive measurements. 



The manuscript is organized as follows: after a descrip- 
tion of the SIAM, we explain the individual features in 
the 4/ density of states and aspects of the numerical NCA 
calculations (Sec. |lT|. In section Hill following section, we 
describe the experimental setup and the sample prepara- 
tion. Section IIVI gives the experimental and theoretical 
results for the different HF compounds, followed by a 
discussion of the quantitative analysis. Finally, in the 
appendix we explain in detail the normalization method 
and the modeling of the spectra. 

II. THEORY OF THE SPECTRAL FEATURES 

The SIAM Hamiltonian for Ce in a metallic host reads, 

ma pma 

+ - V f , ,f , , (1) 

where Hq = X^po- ^p^pa^pa describes the conduction elec- 
tron band, and f}^^ creates an electron with spin cr in a 
4/-orbital with energy Sfm < Ep, m = 1, . . . , 25, S = 
7/2. The 4/-orbitals hybridize with the conduction band 
via the matrix elements Vpm- For later use we introduce 
the effective couplings Tmm' = K'ip^pm'^PT(O) < 

|£/m|, with Apa(uj) the conduction electron spectral func- 
tion and N{lo) — X]p^pcr(^) the conduction electron 
density of states per spin, which in the following will be 
assumed to be flat. The Coulomb repulsion U between 
electrons in any of the local orbitals is large enough to 
suppress any double occupancy of the 4/-shell and may 
be assumed [/ — > oo in the following. For the quantum 
mechanical treatment of spin as well as charge fluctua- 
tions the /-electron operators can then be represented as 
fma ~ 9rna^' whcrc the bosonic and fermionic auxiliary 
operators, b\ gl^ai create the empty and the singly occu- 
pied /-shell (toct), respectively, and obey the constraint, 

Q = Y.9l.9n.. + h%^l . (2) 

ma 

The spectrum of this system has generically six distinct 
features as shown in Fig. [1] (A-F). They can be under- 
stood as follows: At low T the occupation of the lowest 
4/ level is close to unity, n/i < 1, while all other 4/ 
orbitals are essentially empty, n/m « 0, m = 2, . . . ,7. 
Hence, the broad Af^ 4/° ionization peak (A) with a 
fuU width at haff maximum (FWHM) of P w J2m ^im 
corresponds to the lowest single-particle level Sfi. Res- 
onant spin flip scattering of electrons at the Fermi en- 
ergy induces the narrow Kondo resonance (D) of width 
~ fesTft-, shifted hy S ~ sin(7rn/i)A:BTR- above Ep due 
to level repulsion from the single-particle levels e fm (see 
Ref. [2^ for details). The spin-orbit (SO) and the crystal 
field (CF) peaks appear in pairs (B, F) and (C, E), re- 
spectively. They arise from virtual transitions from the 
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FIG. 1: Left: Theoretical 4/ spectrum from NCA calculations based on the SIAM for T= 11 K and the model parameters given 
in the captions of Fig. [5] The hatched region indicates the 'photoemission' region below Ef- The inset shows the near-_BF 
region with Kondo resonance (D) and crystal field features (C,E). B and F are the spin-orbit satellites ( J = 5/2), A is the 
ionization peak. For finite U the two-electron state would appear far right of the displayed energy range at ~ U — |e/|. 
Right: Sketch of the energy level scheme of the SIAM for the conduction band states (c states) and for the impurity (/-states). 



ground state into the (empty) excited SO (F) and CF (E) 
states and vice versa (B and C). The positions of the satel- 
lite peak pairs are, therefore, approximately symmetrical 
about the central Kondo peak (D). However, while the 
features above Ep have significant spectral weight, those 
below Ep appear merely as weak shoulders. This is be- 
cause the transition probabilities carry a detailed balance 
factor w — n^{l ~ n^), where {n-f) is the occupation 
number of the 4/ orbital in the initial (final) state, i.e. w 
is large for the excitations E, F, but small for the transi- 
tions B, C. As n* and n-^ are controlled both by quantum 
and thermal fluctuations the CF and SO satellites are 
T-dependent and are signatures of strong correlations. 

We now analyze the nature of the satellite peaks in 
more detail. A straight-forward Schrieffer- Wolff projec- 
tion onto the subspace of the singly occupied 4/-shell 
yields an s — d exchange model, generalized to multiple 
local levels. 



Hsd - -f^O + ^ >^ fra9lna9raa (3) 
mcr 

pp'rnm' ^i^2^5^i 

cri-t-o"3— o-2-t-(T4 



subject to the constraint ^ with vanishing boson num- 
ber. The effective spin exchange couplings, including 
level renormalizations to second order in Vpm (for non- 



degenerate levels), are obtained asfi^ 



e/1 + E„>i 



(4) 



Due to the anti-commutative operator structure of the 
spin coupling term in the Hamiltonian ([3]) the conduc- 
tion electron-local spin T-matrix acquires in second or- 
der perturbation theory in the Js logarithmic diver- 
gences at the transition energies between the local levels, 

Tmm'{uJ,T) = ~iV(0) JmnJmn' (5) 

n 

{lo + e„, - £„)2 + ttT^ 
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This expression is derived in a straight-forward way us- 
ing the bare auxiliary fermion propagator G^cr{v) — 
{v — \ — em)~^, where the parameter A is taken to in- 
finity to project onto the constrained Hilbert space of 
single occupancy of the 4/-shell (see e.g. Ref. |Tl|). m 
and m! denote which of the 4/-orbitals is occupied in the 
incoming and in the outgoing channel, respectively, uj is 
the energy of the scattering conduction electron, and the 
incoming local particle is assumed to be at the eigenen- 
ergy of the initial orbital, £/m (on-shell). Due to the 
detailed balance factors mentioned above, the divergen- 
cies in Eq. ^ give a significant contribution to the 4/ 
spectral density only, if at least one of the levels m or 



n is the local ground state. The logarithmic energy and 
temperature dependence demonstrates that the SO and 
CF satellite peaks are, in fact, Kondo resonances, i.e. in- 
duced by quantum spin flip processes, however shifted 
by the excitation energies lu — ±(e fm — £ fi) with respect 
to the central Kondo peak near uj = 0. Despite these 
multiple Kondo peaks a single Kondo temperature Tk is 
defined as the crossover scale below which the collective 
spin singlet ground state is formed. Since at tempera- 
tures below the SO and CF splitting energies only the 4/ 
ground state level e fi is significantly occupied, it may be 
estimated as^^ 

Tk « e-i/(2JV(o)j„o) (g) 

and is, hence, given by the width of the central Kondo 
peak near the Fermi energy (D in Fig. [ij. The phonon- 
induced broadening of the Kondo resonance and its satel- 
lites is not considered here, because the coupling of the 
phonons to the spin excitations is small. A detailed anal- 
ysis of the spectral weights and widths of the multiple 
resonances will be given elsewhere, see also Ref. [1^. 

The self-consistent NCA equations are formulated in 
terms of the auxiliary fermion and boson propagators 
Gf, Gb, 

[G^^]mm'{t^) ~ {l^ — ^ ~ £fm)Smm' — '^gcrmm'{l^) (7) 

B-\u) = J.- A-£/„0 -SbM , (8) 

where the auxiliary fermion selfenergy Sg is a matrix in 
orbital space because of the non-conservation of the or- 
bital degree of freedom m, Smm' denotes the Kronecker 
delta and the superscript —1 matrix inversion. In NCA 
the selfenergies Eg , Ef, and the physical 4/ spectral func- 
tion Af are given by, 

^gamm'ii^) = Fw J de [1 ~ /(£)]iV,(e)B(i/ - e) (9) 

)(io) 

(7 7nm' 

A,._,H=/dse-^^^[Ag.w(c. + s)^(e) 

+Ag„,nm'ie)Ab{e-Lj)] (11) 

where all propagators are understood as the re- 
tarded ones, and the spectral functions are Aglv) — 
— ImG(i/)/7r, etc. See Refs.fTlllTslfor the exact projection 
onto the physical Hilbert space of no multiple 4/-shell oc- 
cupation and for an efficient numerical evaluation of the 
NCA equations. 

The NCA is known to correctly include the logarith- 
mic perturbative corrections in the absence of a mag- 
netic fiel d^^i^^ and to reproduce the correct Kondo peak 
width,— as long as the temperature is not too far below 
Tk- It may also be shown to incorporate the above- 
mentioned detailed balance factors and, hence, describes 
well the spectral features discussed above (Fig. [5]). 




energy co [meV] 

FIG. 2; Theoretical 4/ spectrum from NCA calculations 
based on SIAM for different temperatures. The displayed en- 
ergy range covers the Kondo resonance close to zero and the 
crystal field structures at approx. —10 meV and -1-20 meV. 
The inset shows the Kondo resonance at low temperatures, 
i.e. at T = Q.ITk (unitarity limit). 



In order to estimate Tk of the experimental results by 
fitting of the NCA spectra, one calculates the respective 
NCA spectra at T = O.ITk « 1 K, i.e. close to the 
unitarity limit. Tk can then be determined from the 
peak width of this theoretical low-T spectrum (see inset 
ofFig.lH). As already seen in Fig.[I]at low T, i.e. T « Tr-, 
the central KR appears as a narrow line, well separated 
from the CF satellites above and below the Fermi level. 
Towards higher temperatures, the line width of the KR 
increases, while the maximum intensity becomes smaller 
and the distinction of the KR and the CF features is 
successively smeared out. At T Tk the KR disappears 
as separate peak and an enhanced 4/ density of states 
(DOS) near Ep is rather due to the CF excitations. Thus, 
the persistence of an enhanced DOS at the Fermi level 
even at high T in comparison to Tk is naturally explained 
within the SIAM in combination with CF excited states 
and their logarithmically wide extension towards high T. 



III. EXPERIMENTAL SETUP 

The PES experiments have been performed with a 
SCIENTA SES 200 analyzer in combination with a 
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monochromatized GAMMADATA VUV lamp at photon 
energies of hi' = 21.23 eV (He !„) and hiy = 40.8 eV 
(He Hq). The base pressure of the UHV system was 
below 5 X 10~^^ mbar, increasing during the measure- 
ments — due to the He leakage from the discharge lamp 
— to < 1 X 10^^ mbar. The samples could be cooled 
down to approximately T = 4 K on the manipulator in 
the spectrometer chamber. For the presented data anal- 
ysis the calibration of the spectrometer is very crucial, in 
particular the sample (surface) temperature, the energy 
resolution, and the position of the Fermi level. For this 
purpose we have repeatedly performed low-temperature 
reference measurements on polycrystalline Ag and on 
poly-crystalline Nb in the superconducting state. By an 
analysis of these data one can independently determine 
the required parameters. ^^'^^ For the measurements pre- 
sented here the energy resolution of the spectrometer was 
chosen to 5.4 meV as a compromise between intensity — 
in particular when using He 11^ radiation — and energy 
resolution. More about the spectrometer can be found in 
Ref.[27| 



A. Sample preparation 

The five different Ce compounds investigated in this 
work belong to the class of HF compounds which 
are characterized by large specific heat coefficients 70. 
Whereas CeCug displays with 70 = 1.6 JK"^mol~^ 
the largest value^^ among the investigated compounds, 
the other systems range from 70 — 1.1 JK^^mol^^ 
for CeCuzSia^i^ down to 70 = 0.104 JK-^mol-i 
for CeSi2^ (CeRuaSia^^ CeNiaGez^: 70 = 

350 mJK~^mol^^). None of these compounds shows 
any magnetic order down to temperatures of a few 
mK,a>33^i35^i31.38^ In the case of CeCusSia and 
CeNi2Ge2 there is a superconducting phase transition at 
rc = 0.5 and 0.1 K'^a, respectively. 

The single-crystalline CeCug samples were grown by 
Czochralski technique in a high-purity argon atmosphere 
using a tungsten crucible. In contrast to this the poly- 
crystals CeSi2, CeCu2Si2, CeRu2Si2, and CeNi2Ge2 were 
produced by first melting the respective stoichiometric 
ingredients also under high-purity argon atmosphere and 
then tempering in a dynamic vacuum. The purity of the 
used elements amounts to 99.99 % at least. The tem- 
perature and the duration of the temper process varies 
from system to system between 800-1200°C and 50- 
120 hours. All these compounds were characterized as 
single phase compounds by the use of x-ray diffraction 
technique, that also yields the crystal structure and lat- 
tice parameters. Whereas all the ternary compounds 
crystallize in the tetragonal ThCr2Si2-structure^i'^i^ 
with lattice constants a = 4.10 A, 4.20 A, and 4.15 A, 
and c = 9.93 A, 9.80 A, and 9.85 A for CeCuaSia, 
CeRuaSia, and CeNiaGca, respectively, CeSi2 crystalhzes 
in the tetragonal a-ThSia-structure— with a = 4.19 A 
and c = 13.91 A. The crystal structure of CeCug under- 



goes a change from orthorhombic symmetry at tempera- 
tures above 200 K to monoclinic at low temperatures 4^ 
Because the following investigations on this compound 
have all be done in the monoclinic phase, only the lat- 
tice parameters for this structure are given: a = 5.08 A, 
6=10.12 A, c = 8.07 A, and /3 = 91.36°. 

To prepare clean surfaces the samples, cut to a typical 
size of 2 X 3 X 5 mm'^ and equipped with a cleavage post, 
were fractured in situ at low temperatures just before the 
UPS measurement. The resulting surfaces were usually 
coarse grained, even when the sample was single crys- 
talline. In the latter case, the rough surface topology 
ensures an integration over a sufficiently large effective 
fc-space. This is important, because a strong angular de- 
pendence of the 4/ intensity was observed4^i^ Scraping 
the surface with a diamond file did not yield satisfying 
results.«4S' Because of the high surface reactivity of the 
rare-earth compounds the duration of the measurement 
was kept below 12 h. During this period the surface qual- 
ity was repeatedly checked by a measurement of the O 2p 
photoemission intensity at binding energies of about 6 eV 
in the valence band spectra. 



IV. RESULTS AND DISCUSSION 

In the following section we present our results on five 
HF compounds with Kondo temperatures below 100 K, 
namely CeCug {Tk = 5 KiS), CeCu2Si2 (10 K^O), 
CeRuaSia (20 K^), CeNi2Ge2 (30 K^), (41 K^); the 
given Kondo temperatures Tk are determined by inelas- 
tic neutron scattering (INS). Here we restrict our inves- 
tigations to Iow-Tk or 7-Ce like materials with a low 
hybridization strength because of two reasons: First it 
is a well known fact that these 7-like materials only ex- 
hibit a small change in hybridization strength by going 
from the bulk to the surface that is not the case in a-like 
Ce compounds.^— Recently J. W. AUen^^ showed that 
PE spectra of 7-like Ce compounds mostly exhibit bulk 
rather than surface contributions. Second, the position 
of the KR in heavy fermion Ce compounds occurs above 
Ep at (5 w sminn fi)kBTK -^^ Because our normalization 
method (described in detail in Appendix [A]) only allows 
to observe spectral features in the energy range between 
Ep and bhsT above Ep, it is a crucial aspect to find the 
KR in this energy interval at low temperature. By using 
the normalization procedure we are also able to compare 
the PE spectra with NCA spectral functions. To get 
a quantitative comparison, we have to apply this nor- 
malization method to both the PES data and the NCA 
results. By the performance of a conventional x^-fitting 
procedure, we get a suitable parameter set for the de- 
scription of the spectral function of Ce systems as de- 
scribed in Sec. |TT1 
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FIG. 3: Photoemission spectrum of CeCug (left panel: He IIq with hv = 40.8 eV; right panel: He la with hv = 21.2 eV). 
Filled circles represent the raw data {AE = 5.4 meV, T = 16 K), mainly consisting of the tail of the KR. The SO partner 
(J = 7/2) appears at a higher binding energy of ~ 250 meV given in the inset {AE « 15 meV). Open circles represent the 
normalized data using the experimentally broadened FDD, the shaded area marks the unreliable spectral range above 5fcsT. 
The solid line represents the fitted NCA spectral function with e/ = — 1.05 eV, D — 2.8 eV, CF splittings of the J — 5/2 sextet 
AcF = 7.2/13.9 meV, SO splitting Aso = 250 meV, hybridization V^lld meV. 
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TABLE I: Comparison of the Kondo temperatures Tk and the CF energies Acf determined from PES and from other 
experimental methods, namely from (a): INS studies, (6): specific heat measurements, (c): Raman scattering experiments, and 
(d): theoretical considerations based on specific heat measurements. 



A. Low-temperature spectra of CeCug 



The first presented system is the prototype HF com- 
pound CeCug. Because of the small Kondo-temperature 
of only Tk ~5 K and the correspondingly weak 4/ spec- 
tral weight at the Fermi level, there exist only few valence 



band photoemission studies on this compound i^'^'^'^ 
Our high resolution photoemission spectra on CeCug are 
displayed in Fig.[3]measured at r= 16 K with both He la 
(right panel) and He IIq. (left panel) radiation. The insets 
show an extended energy range below Ep over 500 meV 
with the weak tail of the Kondo resonance at Ep and the 
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FIG. 4: Temperature dependence of the experimental (PES) and calculated 4/ spectra of CeCug after the application of the 
normalization procedure. The PES data are taken with He IIq (left panel) and He 1^ radiation (right panel). The used NCA 
model parameters are given in the caption of Fig. [S] The insets show the raw spectra, normalized to the same intensity at a 
binding energy of « 100 meV. 



SO feature at approx. 260 meV corresponding to peak B 
in Fig. [T] Although the photoionization cross section of 
the Ce 4/ states increases about a factor 3.2 from an exci- 
tation energy of hv= 21.23 eV to 40.8 eV^"* both spectra 
clearly reveal the 4/ derived spectral features over the 
conduction band background, although the 4/ intensity 
with He I is significantly reduced. 

In the blowup of the neax- Ep region in Fig. [31 only 
the tail of the KR is discernible, with an increasing in- 
tensity towards Ep. The spectral information above the 
Fermi level — the maximum intensity of the Kondo res- 
onance is anticipated slightly above Ep — is suppressed 
by the Fermi-Dirac distribution (FDD) that steeply goes 
to zero for small temperatures. The spectral information 
can partly be restored when the above mentioned nor- 
malization method (see Appendix |^ is applied to the 
raw He data. As dernonstrated already for He Il-data 
on CeCu2Si2 in Ref. 16, there appears obviously a nar- 
row peak with a full width at half maximum (FWHM) 
of w 6 meV and a maximum at about 3 meV above Ep. 
The peak can be observed in both the He 11^ and in 
the He Iq spectra, but with a peak-to-background inten- 
sity ratio which is about a factor of 2 larger for He II. 
The peak maximum is sufficiently below the upper limit 
of bksT « 7 meV for the applicability of the method 
bksT « 7 meV, where the noise of the data becomes 
enormously magnified due to the exponential decrease of 



the FDD. It should be mentioned that the data for bind- 
ing energies below ikgT (below Ep) remain unchanged 
by the normalization procedure. 

In order to analyze the narrow peak feature further, 
we have performed NCA calculations with a set of model 
parameters that is given in the figure caption of Fig. [3l 
The starting points for the model parameters are taken 
from the photocmission spectra at higher energies (e/, 
A50) or from other experiments, as e.g. inelastic neu- 
tron scattering (Ac_f), and then reasonably modified to 
fit the normalized experimental spectra. The intensity 
ratio between 4/ and conduction band states, which are 
assumed to be constant in energy, is set arbitrarily, the 
Coulomb energy U is infinite. The agreement between 
theory and experiment is striking. The NCA result is 
able to describe exactly the experimental line shape and 
the energetic position of the narrow peak just above Ep 
both in the spectra taken with He 11^ and He radia- 
tion. 



B. Temperature dependence of the CeCus 4/ 
spectra 

The NCA allows to calculate spectra for finite tem- 
peratures. All physical properties — and in particular 
the spectral density of statesi^ — scale with the Kondo 
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FIG. 5: Photoemission spectra of CeCu2Si2 (He 11a, hi/ = 40.8 eV). Left panel: normalization at T = ll K (cf. Fig. |3]); the 
inset shows an extended energy range (A_E « 15 meV) including the SO (J = 7/2) excitation at ~ 270 meV below Ep taken 
with He 11c, radiation (open circles) and with He Iq, (filled circles). Right panel: temperature dependence of the experimental 
(He Ha) and the calculated 4/ spectral function of CeCu2Si2 after the normalization. Model parameters: e/ = —1.57 eV, 
Z) = 4.3 eV, AcF = 32/37 meV, Aso = 270 meV, 1^ = 200 meV. The inset shows the raw data normalized to the same intensity 
at a binding energy of ~ 100 meV. 



temperature Tk- In addition, the recoverable energy 
range (above Ep) increases linearly with T, because the 
thermal broadening of the FDD increases with the sam- 
ple temperature. Therefore we analyzed the near-E'i? 
spectral structures of CeCue for different temperatures. 
Fig. m shows the experimental data plus the NCA result 
at several temperatures, normalized at the FDD at the 
respective temperature (see Appendix \^ . Again we ap- 
plied the normalization to both the He IIq (left panel) 
and the He Hq. (right panel) spectra on CeCug. The in- 
sets give the raw data, normalized at the same intensity 
at Eb = 100 meV. These data are dominated by the in- 
creasing broadening of the FDD and the vanishing tail of 
the Kondo resonance. 

As already seen in Fig. [3] at low T, the KR appears 
in the recovered data at r= 16 K as a narrow line just 
above Ep. Towards higher temperatures the line width of 
the KR increases while the maximum intensity becomes 
smaller. The spectra of systems which clearly resolve 
a crystal field (CF) structure (see below CeCu2Si2 and 
CeNi2Ge2) show an additional smearing of the crystal 
field fine structure. At T = 101 K and above, the narrow 
peak of the KR has disappeared, and is replaced by a 
broad structure with a maximum at w 10 meV above 
Ep, which results in a considerable 4/ density of states 



right at the Fermi level. 

With the same model parameters as for the low- 
temperature spectrum in Fig. [3] we have calculated the 
spectra at the higher experimental temperatures. The 
temperature dependence observed in the experimental 
data is perfectly reproduced by the NCA result which 
is qualitatively described in Fig. [2l As we will discuss 
further below, the existence of the crystal field splittings 
of the 4/ levels has a considerable influence on the dis- 
tribution of spectral weight near the Fermi level. 



C. CeCu2Si2, CeRu2Si2, CeNi2Ge2, CeSia: Low-T 
spectra and temperature dependence 

This becomes more clear in the case of the HF- 
compound CeCu2Si2 with a Kondo temperature of Tr- w 
10 K and a crystal field energy of about 30 meV (see 
references given in Tab. [J). At Eb ~ 20 meV the raw 
low-temperature spectra show a clearly distinct crystal- 
field structure in the 4/5/2 levels<^ A similar CF struc- 
ture could also be observed in CeBg^ and is also visi- 
ble in the CeNi2Ge2 spectra below. In addition to this 
fine structure, the high-temperature data show a large 
shift of the broad peak maximum to approx. 45 meV 
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energy below [meV] energy below [meV] 



FIG. 6: Photoemission spectra of CeRu2Si2 (He Ila, hv = AQ.9, eV). Left panel: normalization at T = ll K (cf. Fig. |3]); the 
inset shows an extended energy range (A_E « 15 meV) including the SO (J = 7/2) excitation at ~ 290 meV below Ep taken 
with He Ha radiation (open circles) and with He Iq, (filled circles). Right panel: temperature dependence of the experimental 
(He Ha) and the calculated 4/ spectral function of CeRu2Si2 after the normalization. Model parameters: e/ = —1.54 eV, 
D = A.2 eV, AcF = 18/33 meV, Aso = 295 meV, 1^ = 170 meV. The inset shows the raw data normalized to the same intensity 
at a binding energy of ~ 100 meV. 



above Ep, even a two- fold fine structure in the high- 
temperature peak can be resolved. Both features, the 
CF satellite below Ep in the low-temperature data and 
the broad split peak at high temperatures are well de- 
scribed by the NCA calculation including the crystal field 
resolved degeneracy of the 4/ levels. Without this CF 
splitting, there is only the spin-orbit splitting on the scale 
of Aso = 250 meV, which does not contribute to the in- 
vestigated energy range, and the observed temperature 
dependence of the spectra can not be described in the 
SIAM quantitatively. The high spectral intensity at Ep 
at temperatures T ^ Tr- is an immediate result of the 
existence of the crystal field structures. 

We performed analogous investigations on the HF com- 
pounds CeRu2Si2, CeNi2Ge2, and CeSi2 which exhibit 
Tr-'s in ascending order from 16 to 40 K. The results 
of our PES investigations are shown in Figs. [BHS] In 
each figure the left panel shows the analysis of the low- 
temperature data and the right panel gives the tempera- 
ture dependence of the weai-Ep spectral features. In the 
following we restrict ourselves to the He IIq data because 
all 4/ states are much more pronounced than at low pho- 
ton energies (He I). Especially in the case of the ternary 
Ce compounds displayed here, the intensity of the 4/ 
spectral features below Ep is reduced dramatically, in 



contrast to CeCug for which the 4/ features could still 
be observed in the He I data. Even by irradiation with 
He Hq the near-Ep spectra of CeRu2Si2, CeNi2Ge2, and 
CeSi2 show only a weak 4/ fine structure. 

A comparison of the He I and He II spectra is given 
in the insets of Figs. O-jS] Taking the difference of the 
two spectra^^ is a useful method to extract the 4/ con- 
tribution from other states near Ep. This can be partic- 
ularly important, when there is a strong modulation of 
the non-4/ background intensity, which is e.g. the case 
for CeRu2Si2 and CeNi2Ge2^^~^^. Here one has the over- 
lap with Ru Ad and Ni Zd bands, respectively, that have a 
comparatively high photoemission cross-section for pho- 
ton energies in the UV rangCfSi In the case of CeCug, 
CeCu2Si2, and CeSi2 the contribution of non 4/ states 
in the interesting energy range is rather small and — 
although not negligible — without a fine structure on 
the scale of a few tens of meV. Another method to ex- 
tract the 4/ spectra from the net photoemission signal 
is resonant photoemission spectroscopy (RESPES) at the 
Ce Sd or 4(i absorption edges Although this method 
is limited in energy resolution and therefore not suitable 
for investigations of the fine structure on the meV scale, 
the relative intensity and energetic position of the ion- 
ization peak (4/°) can be investigated in detail. For all 
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energy below [meV] energy below [meV] 

FIG. 7: Photoemission spectra of CeNi2Ge2 (He IIq, hv — 40.8 eV). Left panel: normalization at T = 11 K; the inset shows 
an extended energy range (A_B«15 meV) including the SO {J — 7/2) excitation at ~ 300 meV below Ef taken with He IIq 
radiation (open circles) and with He I^ (filled circles). Right panel: temperature dependence of the experimental (He IIo,) 
and the calculated 4/ spectral function of CeNi2Ge2 after the normalization. Model parameters: e/ = — 1.43 eV, D — 3.9 eV, 
AcF = 26/39 meV, Aso = 275 meV, 1^ = 209 meV. The inset shows the raw data normalized to the same intensity at a binding 
energy of ~ 100 meV. 



the systems measured in this work we find values for the 
/^-energy (partly from hteratur o^^i^^'^° ) that amounts to 
1.7 eV, 2.6 eV, 2.2 eV, 2.2 eV, and 2.35 eV for CeCug, 
CeCu2Si2, CeRu2Si2, CeNi2Ge2, and CeSi2, respectively, 
similar to the values from non- resonant PES.-^ As de- 
scribed above, we used these values as benchmarks for 
the e / values in the parameter sets of the NCA calcula- 
tions (see figure captions Figs.[5H5]). 



D. Quantitative spectral analysis: Fitting with 
NCA 

As described in Appendix [X] we iteratively fit the nor- 
malized photoemission data for a quantitative data anal- 
ysis. This analysis yields one characteristic model param- 
eter set for each compound, which is used for the calcula- 
tion of the spectra at all investigated temperatures. The 
resulting parameters are given in the captions of Figs. [SHE] 

In addition, each parameter set can be assigned to a 
certain Kondo temperature Tr-, which is the Kondo tem- 
perature of the respective system. To get the Tk values 
we calculated the NCA spectra for each parameter set in 
the unitarity limit at T = 0.1 x Tk- The linewidth of the 
Kondo resonance in this low-temperature limit defines 



immediately the Kondo temperature by FWHM~ /cbTr- 
due to the scaling properties of the SIAM. The resulting 
Tff's are given in Tab. [T] together with the crystal field 
energies A^f and the Kondo temperatures determined 
from other experimental methods, i.e. INS and transport 
measurements. One should note that these reference bulk 
values are spread over a wide range, indicating some sys- 
tematic uncertainties in the determination. 

For all compounds we find an excellent quantitative 
agreement of our Tk values with the Kondo temperatures 
determined by INS studiesi^'^i^^'^^'^^-^'' In addition, we 
find a good coincidence of the crystal field splittings for 
CeCufi, CeCu2Si2 and CeSi2 with the INS values from 
Refs. IsHlll. In the case of CeRu2Si2 and CeNi2Ge2 the 
CF energy separations could not be determined exper- 
imentally upto now, mainly because of uncertainties in 
the assignment of the CF level schemei^i^ Therefore in 
the case of CeRu2Si2 we compare our values with theo- 
retical results proposed by Zwicknagl et al.r^ based on 
specific heat measurements^ and again we find an excel- 
lent agreement. A recent study on CeNi2Ge2 claims the 
first observation of CF structures in INS experiments on 
this compound, giving values that correspond rather well 
with our results 1^ 

In summary, the comparison shows that all of our re- 
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energy below [meV] energy below [meV] 

FIG. 8: Photoemission spectra of CeSi2 (He 11^, hv = eV). Left panel: normalization at r = 11 K; the inset shows an 
extended energy range (A_E«15 meV) including the SO {J = 7/2) excitation at « 280 meV taken with He Ha radiation (open 
circles) and with He lo, (filled circles). Right panel: temperature dependence of the experimental (He Hq,) and the calculated 
4/ spectral function of CeSi2 after the normalization. Model parameters: e/ = —1.35 eV, D = 3.7 eV, Acf = 25/48 meV, 
Aso = 270 meV, V = 203 meV. The inset shows the raw data normalized to the same intensity at a binding energy of 
« 100 meV. 



suits are in reasonable agreement with the values from 
the other methods. This is surprising because of two pos- 
sible problems in our method: 1.) In contrast to trans- 
port and INS measurements is PES a surface sensitive 
technique and the surface properties of a rare-earth com- 
pound are not necessarily identical with the bulk prop- 
erties, in particular the Kondo temperatures might be 
completely different. 2.) The SI AM, on which the NCA 
calculations of the spectra are based, contains several 
significant simplifications that might influence a quan- 
titative description of width and position of the Kondo 
resonance. 

Considering the first point, one has to mention that 
there is an ongoing debate in how far PES with excita- 
tion energies in the VUV-range will be able to reflect the 
hulk properties of Ce compounds. At these photon ener- 
gies (and the respective photoelectron kinetic energies in 
the range of 40 eV for He II) the information depth is at a 
minimum of a few lattice constants.— Therefore, one per- 
forms additional photoemission measurements at higher 
photon energies to increase the information depth. PES 
with photon energies in the soft x-ray range, e.g. in reso- 
nance with the Ce 3d-4/ absorption edge at hv « 880 eV, 
is regarded as a sufficiently bulk sensitive method with 
an information depth of about three times larger than 



for He IIi^ On this basis Sekiyama et al^^i^^^ have 
performed 3d-4/ RESPES studies on CeRu2Si2 with a 
noticeably high energy resolution (Ai?== 120 meV). They 
compared the experimental spectra with results from 
NCA calculations, including CF splittings of 51.6 meV 
and 68.8 meV,— and found a good coincidence between 
the RESPES data and the NCA spectrum. The signif- 
icant difference between these crystal field energies and 
our values prompt us to investigate this point further. 
We found the following results: 

(i) If we perform NCA calculations with Sekiyama's 
parameter set at T = 20 K we find clearly separate 
CF peaks below and above Ep. As seen in Fig. [5] 
we positively cannot observe such structures below 
Ep in our high-resolution low-temperature spectra. 
In addition, these CF energies would be too large 
to reproduce the observed temperature dependence 
in Fig. [SI right panel. 

(u) If wc convolute the NCA spectra at T = 20 K 
with our parameter set (see captions Fig. [5]) with 
a Gaussian to describe an energy resolution of 
AE = 120 meV, we can reproduce Sekiyama's 3d- 
4/ RESPES spectrum with the same accuracy. 

Therefore we conclude, that the RESPES results on 



12 




energy above [eV] energy above E^, [eV] T [K] 

FIG. 9: Left panels: NCA spectra for CeSi2 (parameters given in the captions of Fig. (a) Raw and (b) convoluted with a 
Gaussian of FWHM= 600 meV. Panel (c) shows the (normalized) temperature dependence of the /^-intensity as calculated from 
the convoluted NCA-spectra in panel (b), in comparison with the results of IPES-studies.— Panel (d) shows the temperature 
dependence of the 4/ occupation n/ as resulting directly from the NCA calculations. These values are compared with the ones 
extracted from XAS measurements at the Ce L3 edge.— 



CeRu2Si2 do not contradict our results and, further more, 
do not indicate a quantitative deviation of the crystal 
field energies extracted from surface sensitive from the 
bulk values. 

There is another possible consistence check of our 
analysis, using the comparison with temperature de- 
pendent (resonance) inverse photoemission spectroscopy 
(RIPES), which immediately has access to the spectra 
above the Fermi level where for Ce-systems the main 
spectral 4/ weight of the Kondo resonance and its satel- 
lites appears. Again the energy resolution of this method 
is not sufficient to resolve the fine-structure of the 4/^ 
features, namely the Kondo resonance and its crystal 
field and spin-orbit satellites (J=7/2). A second broad 
peak appears around U — |e/| and corresponds to the 
two-electron final state p. However, the NCA spec- 
tra, which we get from our fit of the PES data, contain 
also the spectral information above the Fermi level. For 
the following comparison we have chosen CeSi2 because 
there exist detailed RIPES^^ and x-ray absorption (XAS) 
investigations^^ on the temperature dependence of the 4/ 
occupation number Uf, which can be determined from 
the intensity ratio of the and features. 

To get this temperature dependence from our results 
we model the inverse photoemission spectra by a convolu- 
tion of the NCA spectrum (with parameters from fitting 



the PES data) above the Fermi level (i.e. s{E,T) x [1 - 
f{E,T)] at the respective temperature, cf. Appendix E]) 
with a Gaussian describing the finite energy resolution of 
the RIPES experiment. The maximum position and the 
line width of the result — including Kondo resonance, CF 
satellites and spin-orbit {J = 7/2) partner — are identi- 
cal to the experimental values within the experimental 
errors?^ The left two panels of Fig. [5] shows (a) the un- 
broadened NCA spectra a three different temperatures 
and (b) the result of the convolution (Ai? = 600 meV). 
From this modeled RIPES spectrum we take the tem- 
perature dependent integrated intensity for a comparison 
with the experimental temperature dependence. 

Panel (c) of Fig. [5] gives the temperature dependence 
of the relative integrated RIPES intensities I{f^)/I{p) 
of CeSi2 published in Ref. [zl (open circles) . This inten- 
sity ratio, proportional to 1 — n/, clearly increases at low 
temperatures and the curve even suggests the existence 
of the plateaus well above and below Tk « 30 K as ex- 
pected from the theory. We compare this behavior with 
the integrated intensity from Fig.[5](b), scaled linearly to 
match the RIPES value at the minimum temperature of 
15 K (there is no peak since C/// = 00 in our model) . A 
clear coincidence of the two results is obvious, in partic- 
ular if one takes the error bars of the RIPES experiment 
into account. However, the S-shape of the experimental 
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temperature dependence is less pronounced. 

From an analysis of the NCA spectra one can also 
immediately determine the 4/ occupation number Uf. 
Experimentally, this 4/ occupation can be measured by 
x-ray absorption spectroscopy (XAS), e.g. at the Ce L3 
edge at hi' « 5720 eV. Although this method is truly bulk 
sensitive, one has to know that the final state consists of 
an additional electron in the 5c? valence states, which 
can lead to a slightly modified (final state) 4/ occupa- 
tion sampled by this method. However, Fig. [5] (d) gives 
the temperature dependence of n/ as extracted from an 
XAS experiment at the Ce L3 edge of CeSi2'?® Although 
the shape of both curves is very similar, the XAS values 
are typically 0.01-0.03 smaller than the NCA numbers. 
The authors in Ref. [tI explain this difference by surface 
effects in PES, but this a) can be ruled out by our previ- 
ous observations, and b) is usually exceeded by the final 
state effect in XAS.^^'*° However, the qualitative agree- 
ment between XAS and PES/NCA result is reasonable 
and indicates that our data analysis is reliable. 

Considering the surface effects in the special case of 
CeSi2, which has the highest Tr- of all Ce systems in- 
vestigated in this work, one should note that the pub- 
lished Kondo temperatures, including the values from 
spectroscopy methods, differ significantly in a range of 
« 22 K to 140 K-^L^LI^,^L^ the given Kondo tem- 
peratures of the other compounds have a comparatively 
small scatter. This leads to the impression that CeSi2 
can be regarded as borderline system between 7-Ce- and 
a-Ce like materials, in which the difference between the 
hybridization strengths at the surface (or surface near 
region) and the bulk becomes evident. 



V. CONCLUSIONS 

In the present paper we have demonstrated that high- 
resolution photocmission allows a detailed and quantita- 
tive investigation of the 4/ spectral features of low Tk 
cerium compounds close to the Fermi level. The Kondo 
resonance, which for Ce-systems has a maximum above 
Ep, could be restored for all investigated 7-Ce like HE 
systems by application of a well known normalization 
procedure. The resulting 4/ spectra can be iteratively 
fitted by NCA calculations based on the SIAM, using in- 
dividual model parameter sets including the spin-orbit 
and, in particular, the crystal field one-electron ener- 
gies e/. Our results demonstrate that the consideration 
of the crystal-field splitting and the corresponding fine 
structure in the 4/ spectra is of high importance for the 
consistency of the photocmission results with the ther- 
modynamic properties. Surface effects can be ruled out 
for the investigated 7-Ce like systems. In spite of the 
known limitations of the model, our data analysis yields 
Kondo temperatures and crystal field energies that are 
in surprisingly good agreement with values from other 
experimental methods published in the literature. For 
the investigated temperature range, the description by a 



local model — namely the SIAM — is obviously suffi- 
cient to describe the photocmission data with an energy 
resolution of the order of the Kondo temperature. 
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APPENDIX A: NORMALIZATION PROCEDURE 

On the example of the 3d bands of Ni(lll), Greber— 
has shown that a careful analysis of photocmission data 
close to the Fermi level allows an extended access to the 
thermally excited spectral features up to bksT above Ep. 
Two conditions for this method have to be fulfilled: 1.) a 
satellite free and highly intense photon source must pro- 
vide low noise data, and 2.) the total experimental en- 
ergy resolution AE must be at least comparable to the 
thermal broadening AksT of the Fermi edge; we meet 
these two conditions by using the experimental setup de- 
scribed above. The normalization method has been es- 
tablished for the qualitative spectroscopic investigation 
of the neai-Ep range for various compounds, both for 
angular integrated and angular resolved data (see e.g. 
Refs. ^0,84,85,86) 

The experimental photocmission spectrum / can be 
described as the density of electronic states (or spec- 
tral function) s(E, T) multiplied by the temperature de- 
pendent Fermi-Dirac distribution (FDD) f{E,T). To 
consider the finite energy resolution, the product must 
be convoluted by the spectrometer function g{E, AE), 
which is usually approximated by a Gaussian with the 
full width at half maximum (FWHM) AE. Note that 
in our case the Gaussian describes nearly perfectly all 
reference spectr a^^'^^ , independent from the individual 
spectrometer settings. Thus one gets the photoelectron 
spectrum 

I{E,T) = [siE,T)-f{E,T)]®giE,AE), (Al) 

plus statistical noise, which we neglect in the following 
discussion (see Ref. |8j for more). In principle, this rela- 
tion is also valid for angular resolved spectra, but here 
the finite angular resolution must be taken into account 
too. 

The normalization procedure is realized by dividing 
the measured spectrum I{E, T) by a Gaussian-broadened 
FDD: 

[s{E,T)-f{E,T)](E>g{E,AE) 

f{E,T)®g{E,AE) ' ^ > 
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FIG. 10: Application of the normalization procedure of 
eq. (IA2p at two different temperatures. Top panel: T — 
108 K» AE/kB, lower panel: 9 K< AE/ks- The energy res- 
olution was A_B = 5.4 meV. The grey curves give the modeled 
photoemission spectra, including FDD and Gaussian broad- 
ening, the black solid lines represent the (broadened) intrinsic 
spectral functions. 



Because the convolution is not distributive (i.e. [s ■ f]® 
g ^ [s® g\-[f ® g])^ the result of eq. I|A2|) is not equal to 
the broadened spectral function s{E,T) (g) g{AE,E). In 
other words, the normalization procedure does not a pri- 
ori reveal the intrinsic spectral function. As long as the 
energy broadening AE is small compared to the thermal 



broadening AksT, the deviation between the normaliza- 
tion result and s ^ g is small and the result gives a useful 
qualitative information about the structures in the en- 
ergy range up to bhsT above Ep. 

However, the normalization can be used for even a 
quantitative analysis of the intrinsic spectral function 
when an iterative procedure is applied.— 

The influence of this effect can be demonstrated by 
its application to model spectral functions (similar to 
the theoretical NCA spectra with realistic parameters) 
at different temperatures. The top panel of Fig. [TO] 
gives the broadened intrinsic spectrum s'^^^ ® g with 
AE — 5.4 meV (black solid line) and the modeled pho- 
toelectron spectrum \^s^'~'^ ■ /] ® according to eq. (|A2[) 
(grey line). We divide the latter hy f ®g with T=108 K 
and get the dotted curve, that differs only little from the 
desired result s^'^^(E)g over the investigated temperature 
range of approximately up to ^ksT. 

In this high-temperature case one can replace — as also 
stated in Ref. Issl — the denominator in eq. (|A2[) . i.e. the 
broadened FDD f (g) g, hy a. bare FDD f{E,Teff) with 
an effective temperature of Te// = y^T^ -|- (AE/ikB)^- 
As shown in the top panel of Fig. [TO] this simplified nor- 
malization is in good agreement with s^'^"^ ® g. 

For low temperatures T < AE/ks, the difference be- 
tween the normalized spectrum — with convoluted or 
effective FDD — and the intrinsic spectrum becomes sig- 
nificant (see lower panel of Fig. [TO)) . Both the resulting 
width and the maximum position deviate considerably 
from the right numbers. Surprisingly, the shape of the 
the normalization with / is closer to the intrinsic spec- 
trum s g in the present case than the one with f ® g. 
As described above, the difference between the individ- 
ual curves decreases with increasing temperature, finally 
matching when T >> AE'. 

However, at low temperatures the spectral shape of the 
normalized spectrum must not be used to determine line 
width and position directly. Instead one has to fit the 
spectra iteratively by comparing (e.g. by calculating y^) 
the normalized experimental data with a theoretical spec- 
trum treated the same way. This means one has to model 
the experimental spectrum according to eq. lAU and use 
the same normalization procedure for this modeled spec- 
trum and for the experimental data. If the curves match, 
one can get the position and the line width from the bare 
s^'^"^{E,T) included in the modeled function. 
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